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The measurement of the proton structure function at HERA is often seen as a hint for the 
observation of saturation in high-energy QCD e.g. through the observation of geometric scaling. 
Accordingly, the dipole picture provides a powerful framework in which the QCD-based saturation 
models can be confronted to the data. In this paper, we give a parametrisation of proton structure 
function which is directly constrained by the dynamics of QCD in its high-energy limit and fully 
includes the heavy quark effects. We obtain a good agreement with the available data. Furthermore, 
to the contrary of various models in the literature, we do not observe a significant decrease of the 
saturation momentum due to the heavy quark inclusion. 



I. INTRODUCTION 



The study of the high-energy limit of perturbative QCD has lead to strong predictions for the scattering amplitudes. 
One of the most important result is the property of geometric scaling which is a consequence of saturation [l], 0, EH, 0] 
which extends into the dilute regime where the amplitude is far from the unitarity limit. We shall recall later in this 
paper how this general property can be proven in perturbative QCD from the Balitsky-Kovchegov (BK) equation Q 
or from the Colour Glass Condensate (CGC) formalism 

At small x, the confrontation of those predictions with the experimental measurements of the proton structure 
function at HERA can be achieved within the framework of the dipole model. In the dipole frame, the virtual photon 
fluctuates into a qq pair of flavour / and size r which then interacts with the proton: 



(Q 2 ,x) =V / d 2 r \ dz\^ { L f) T (r,z;Q 2 )\ 2 CT dip (r,x). (1) 
j J Jo 



The photon wavefunction can safely be computed in QED and is found to be 



|v^(r,z;Q 2 )| 2 = e 2 ^AQ 2 z 2 (l-z) 2 K 2 (rQ f ), 

\^ T (r,z;Q 2 )\ 2 = e)^{[zZ + (l-zf]Q)Kl{rQ f ) + m )Kl{rQ f )}, 
where Q 2 = z(l — z)Q 2 + to 2 . The proton structure function, obtained through 

F 2 (Q 2 ,x) = [alP(Q 2 ,x)+a^(Q 2 ,x) 

is thus expressed in terms of the dipole-target cross-section Odip = 2nR 2 T(r, x) where T is the dipole-target scattering 
amplitude as entering the QCD evolution equations. R p is often referred to as the radius of the proton and is to be 
taken as a free normalisation parameter of our model. We are therefore left with the parametrisation of T(r; Y), with 
Y = log(l/a;) called the rapidity. 

This kind of approach is not new in itself. Different approaches to parametrise the dipole-proton scattering am- 
plitude has already been proven successful. One can cite e.g. the pioneering work of Golec-Biernat and Wusthoff 
[tJ which can be improved by adding the collinear DGLAP effects Jl, Q. Those approaches, directly based on the 
gluon distribution function xg(x,Q 2 ), give good description of the data, even including the contributions from heavy 
quarks ■ From a different point of view, some other successful approaches [llj use Regge parametrisations for the 
dipole-proton scattering amplitude. 

Beside those various approaches, perturbative QCD provides definite predictions for T, including its approach to 
saturation. This is perfectly suited for this kind of problem. Indeed, the factorisation formula {J) underlying the 
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dipole picture is valid only in the high-energy (small- a;) limit. Hence, using prediction from QCD at high energy to 
parametrise the dipole-proton scattering amplitude appears as a natural way to proceed. Those predictions that we 
shall recall in is paper have been successfully gathered into a parametrisation for (T<n p (r, x) and tested against the 
HERA data fl2l |. so far including the contributions from light quarks only. All those approach suggest a saturation 
scale Q s , the energy-dependent momentum scale below which the amplitude is saturated, of order 1 GeV 2 for x ~ 
1(T 4 - 1(T 5 at HERA. 

However, it is rather well-known that, once including the heavy quarks (mainly the charm), all approaches observe 
a decrease of the saturation scale. Given those a 'priori important effects of the heavy quarks on the saturation, it 
appears important to reconsider the predictions of the QCD at high energy to include those contributions. In this 
paper, we will obtain two results: first that it is possible to accommodate the predictions from [l2| with heavy-quark 
contributions and, second, that this does not lead to a decrease of the saturation scale. 

In the next section, we shall recall how to build the dipole-target scattering amplitude T from the equations of 
QCD at saturation. We also discuss the idea allowing for the inclusion of heavy quarks. We shall then present the fit 
to the F2 HERA data in itself, including data selection and parameter adjustment. As a conclusion, we shall finally 
discuss our results w.r.t. other models as well as with predictions from NLO BFKL [l3l H4L Il5j]. 



II. QCD PREDICTIONS FOR THE AMPLITUDE AND HEAVY QUARKS 

Since our approach is mainly based on the QCD fit introduced by Iancu, Itakura and Munier (IIM) [T3] , we start by 
a presentation of that model. The formula we use to parametrise the dipole-proton scattering amplitude is obtained 
from our knowledge of the solutions of the BK equation which captures the main ingredients of the high-energy physics 
with saturation effects. The exact solution to that equation is not known but its asymptotic behaviours, for large and 
small dipole sizes has been studied in details. 

In the last years, it has been shown Q that the BK equation lies in the same class of universality than the Fisher- 
Kolmogorov-Petrovsky-Piscunov (F-KPP) equation. The latter has been extensively studied in statistical physics 
over the past seventy years and it is well-known that its solutions can be written in terms of travelling waves. In 
the language of the QCD variables we have used so far this means that if one looks at the rapidity evolution of the 
amplitude T(r]Y) (seen here as a function of r), the amplitude "front" moves towards smaller values of r without 
changing its shape. The "position" of the wavefront is then naturally associated with Q~ , the inverse of the saturation 
scale. At asymptotic rapidities, the amplitude T, initially a function of r and Y independently becomes a function 
of the single variable rQ s (Y). This very important consequence of saturation is the geometric scaling [l7l . Il8l . [l9j 
property which physically means that the physics remains unchanged when one moves parallel to the saturation line. 

This travelling-wave analysis provides two fundamental pieces of information: first, the saturation scale increases 
exponentially with the rapidity 1 : Q 2 S {Y) oc exp(AK). Then, the amplitude is known in the small-r region: 



T(r; Y) oc exp 



/ \ {P-Psf 

-Ic{p-Ps) - 



2a X '<Y 



(2) 



with p = log(4/r 2 ) and p s = log(Q 2 ). The critical slope j c as well as the parameter A = a>x' c and x" are determined 
from the linear BFKL kernel only. This is an important property: though the scattering amplitude fully satisfies the 
unitarity constraints and is sensitive to saturation effects, the parameters which describes it do not depend upon the 
details of how saturation is encoded. We will discuss the value for those parameters later in this section. The fact that 
they do not depend on the details of the saturation mechanism is another interesting feature of this approach. As a 
related comment, one can also obtain [3, HJ the result ([2|) by looking at the BFKL equation with a boundary condition 
at the saturation line: performing a saddle-point approximation gives ((2]) as well as the exponential behaviour of the 
saturation momentum. Note finally that the Gaussian part of the exponential in J5]| violates geometric scaling as it 
introduces an explicit dependence in Y. This term however becomes less and less important as rapidity increases. It 
controls how geometric scaling is approached and allows to deduce that geometric scaling is valid within a window 
P ~ Ps ^ ^/2ax' c 'y. This is an important point that saturation effects are relevant up to large scales above the 
saturation momentum i.e. in the dilute domain. 

The amplitude in the saturated domain is also obtained from the BK equation (2p| (it can also be obtained from 
the Colour Glass Condensate formalism [2l|). Putting it together with equation ^ we reach the final expression for 
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our dipole-proton scattering amplitude: 



T exp 



if rQ s < 2, 



T(r;Y) = { I 2 ^ Y J (3) 

{ 1 - exp [-a(p - p s - b) 2 ] if rQ s > 2, 

where the parameters a and b are fixed so as to ensure that T and its derivative are continuous at rQ s = 2. Going 
from © to (|3]), we have used ax" = (Xc/x' c )^ Y = K ^ Y witn K = Xc/x'c- 

Let us now discuss the parameters in this QCD-saturation-based model. First, the saturation scale can be written 



(Y) = [f) GeV 2 . (4) 

This leaves two free parameters: Xq which is related to the value of the saturation scale at zero rapidity and corresponds 
to the value of x at which Q s = 1 GeV, and A which controls the rapidity evolution of the saturation momentum. 
While leading-order (LO) BFKL predicts A = a\' c ~ 0.9, an analysis of the next-to-leading order (NLO) BFKL 3] 
gives A ~ 0.3, a value which is in much better agreement with the phenomenological analysis. 

Concerning the parameters in the amplitude itself, we shall fix the matching amplitude To. As in [l2|, a default 
value of 0.7 gives good results and variations around that value leads only to small differences, so we will adopt 
T = 0.7. The value of k will be set from the LO BFKL kernel which gives k « 9.9. Note that the NLO BFKL 
predictions, though a little bit smaller, remain of the same order. 

The value of the critical slope 7 C is a fundamental issue of this paper. In the original work where only the 
light quarks were considered, the authors has fixed it to the value obtained from the LO BFKL kernel (j c ps 0.6275). 
However, as we shall see in section IIII1 when including the heavy quarks, keeping that value leads to a dramatic 
decrease of the saturation momentum together as well as to a poor \ 2 for the fit. A key issue of the present work is 
to show that, allowing that parameter to vary, we recover a similar saturation scale and a good fit. In addition, we 
shall see that the value for j c coming out of the fit is rather close to what we expect from NLO BFKL (-f c > 0.7). 

The last parameter entering the dipole-proton cross-section (Tdip(r,x) is the radius of the proton R p which fixes 
the normalisation w.r.t. the dipole-proton amplitude T(r,x). We are thus left with 4 parameters: R p , xq and A 
which were already present in [12|], and the new one: j c . In the massless case, it was found, for To = 0.7, A = 0.253, 
x = 2.67 10~ 5 and Rp = 0.641 fm for 7c fixed to 0.6275. 

Last but not least, we have to specify our heavy-quark prescriptions. While in [l2j the sum in ([TJ) was restricted 
to the three light flavours only, we now want to consider the effect of massive charm and bottom quarks. Those 
are thus introduced in the sum over all flavours. The quark masses, entering the photon wavefunction are fixed to 
mu,d,s = 140 MeV, m c — 1.4 GeV and nib = 4.5 GeV and we have used the modified Bjorken variable x(l + 4m 2 / Q 2 ) 
in the contribution of the heavy quarks. Note that the contribution of the charm and bottom quark to |T]) directly 
give the charm and bottom structure functions. 



III. FIT TO THE HERA DATA 



Now that we have fully described the saturation-based model we are using to describe the DIS structure function, we 
will test its validity by fitting its free parameters to the experimental measurements of the proton structure function 
F%. In this section we present the details of the fit and discuss the results. 

We first have to specify which dataset we are working with. Following the recent analysis, we shall use the last 
HERA data i.e. the last ZEUS[23j and HI [13] measurements of F2. We include a 5% renormalisation uncertainty on 
the HI data to account for a normalisation discrepancy between ZEUS and HI. Note that the analysis in [12] only 
takes into account the ZEUS data. We shall come back on this point later, when we turn to the discussion of our 
results. 

Our approach is focused on the high-energy behaviour of DIS: the dipole-model factorisation [T] is only valid at 
sufficiently small x and, accordingly, the dipole-proton amplitude is build from the high-energy QCD equations. 
Hence, we shall limit ourselves to x < 0.01, the usual cut in those approaches. Furthermore, our approach, based on 
the QCD equations describing saturation, does not takes into account the DGLAP corrections beyond the double- 
logarithmic approximation. We shall restrict our analysis to Q 2 < 150 GeV 2 . Note that the former analysis with 
massless quarks [l2j is more conservative and uses Q 2 < 45 GeV 2 . However, in order to check the validity of the 
results discussed hereafter, we have tested both cuts and observed that the \ 2 and the parameters were not significantly 
changing when increasing the Q 2 cut. 

We have fitted the free parameters of the model discussed in Section [H] to the 281 data contained in our dataset. In 
order to grasp the effect of the correct treatment of the heavy-quark masses, we performed the fit with and without 
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7c 


A 


xo (10" 4 ) 


R P (GeV- 1 ) 


X 2 /n.o.p. 


light quarks only 


[12] (ZEUS only) 


0.6275 


0.253 


0.267 


3.25 


- 


7 C fixed 


0.6275 


0.2574 ± 0.0037 


0.2750 ± 0.0240 


3.241 ± 0.018 


0.959 


7c free 


0.6194 ± 0.0091 


0.2545 ± 0.0051 


0.2131 ± 0.0651 


3.277 ± 0.044 


0.956 


light+heavy quarks 


7c fixed 


0.6275 


0.1800 ± 0.0026 


0.0028 ± 0.0003 


3.819 ± 0.017 


1.116 


7c free 


0.7376 ± 0.0094 


0.2197 ±0.0042 


0.1632 ±0.0471 


3.344 ± 0.041 


0.900 



TABLE I: The table gives the parameters and \ 2 P er point obtained from the fit. The results are shown with or without the 
heavy quark contribution and with 7 C free or fixed. The last line of this table is the main result of this paper. 




FIG. 1: Result of our fit to the proton structure function. The left plot shows the low Q 2 bins (Q 2 < 1 GeV 2 ) while larger Q 2 
are given on the right plot. The Q 2 value corresponding to each curve is given (in GeV 2 ). For clarity, successive curves, from 
bottom to top, have been rescaled by powers of 1.5. 

including the charm and bottom contribution to F2. In addition, for both cases we give the result of the fit for the 
critical slope 7 C fixed to its LO value or considered as a free parameter. The resulting parameters and \ 2 are presented 
in table [J where we have also added the initial parameters from [l2j (with ZEUS data only and Q 2 < 45 GeV 2 ) for 
better comparison. 

Those results deserve some comments: 

• Concerning the re-analysis of the fit without heavy quarks (second line of table [J), we see that with the addition 
of the HI data and the extension of the Q 2 domain, the parameters remain similar and the x 2 is good. 

• If one allows j c to freely vary in the massless case (third line of table [TJ) , again, the fit naturally converges to a 
minimum which is close to the initial one without improving significantly the x 2 . The LO choice 7 C = 0.6275 is 
even compatible with the error bars of what we obtain when we fit it. 

• Once the heavy quarks are taken into account, the situation changes drastically. If the critical slope is fixed to 
its LO value, the situation becomes dramatic. Indeed, not only the quality of the fit is getting worse, but also 
the saturation scale is going down by two orders of magnitude (the exponent A is also decreasing significantly) . 
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limit of geometric 
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FIG. 2: For both the fit in this paper (light+heavy flavours, solid lines) and the fit of [12j (light flavours only, dashed lines), 
we plot two information: the saturation scale Qf (x) and the extension of the geometric scaling window log(Q 2 /Qj) < V2kAY. 
We see that the inclusion of the heavy quarks does not lead to a strong decrease both those scales. The points on the plot 
represent the (x,Q 2 ) position of the HERA measurements. 

• If as anticipated in Section [TT| we allow the critical slope j c to vary, the fit (last line of table fl} converges back 
to a good description (even a better \ 2 than the corresponding massless fit). The parameters describing the 
saturation scale are also rather close to those obtained in the massless case. 

• The value of the critical slope we obtain from the fit seems much larger than the LO result used previously. 
However, if one extract that value from various renormalisation-group-improved NLO BFKL kernels [l5[ one 
get a value of 7 C slightly larger than 0.7. This is again in good agreement with the value to which the new fit 
naturally converges. 

• To test the dependence upon the fixed parameters of our model, we have performed various fits varying those 
parameters around their default value. We have thus changed the Q 2 cut from 150 to 45 GeV 2 , included only 
the heavy charm (without bottom), varied the masses of the heavy quarks and varied the matching point To. 
For all those variations, both the quality of the fit and the values of the parameters remained similar, which 
enforces the robustness of the present parametrisation. 

The last line of table [H giving a description of the saturation effects in DIS including the heavy quark effects, has 
to be considered as the main result of this paper. The corresponding description of the proton structure function is 
plotted in figure Q] 

One of the most interesting point of this parametrisation is that the saturation scale obtained with heavy quarks 
included is very similar to the one obtained with light quarks only. This is better seen in figure [5] where we have 
plot the saturation scale (lower curves) as well as the limit of the geometric scaling window both with and without 
heavy quarks contributions. One clearly see that the addition of the heavy quark contribution only slightly reduces 
the effect of saturation. This is an important result as previous models including heavy quark effects all report a 
decrease of the saturation momentum by (roughly) a factor of 2 once those heavy quarks contributions are included. 
From figure [5] it also appears that a large number of data (all data from small Q 2 up to the limit of the geometric 
scaling window) are sensitive to saturation effects. 

Finally, we can compare the predictions of our parametrisation with the HERA measurements [24, 25] of the charm 
and bottom structure functions. Those are naturally obtained from our formalism by taking the charm or bottom 
contribution to the photon-proton cross-section JT]) . The prediction for our model are plotted in figures [3] and [4] for 
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FIG. 3: Predictions of our fit for the charm structure function. 
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FIG. 4: Predictions of our fit for the bottom structure function. 



the charm and bottom structure functions respectively. In both cases, we observe a good agreement with the data. 
Similarly, by taking the contribution coming from the longitudinal part of the wavefunction in {!]), we can obtain 
predictions for the longitudinal structure function. Our result is shown in figure [5] together with the HI measurements 
|22| . Again, the present parametrisation gives a good description of the data. 



IV. DISCUSSION AND CONCLUSIONS 



In this paper, we have shown that it was possible to accommodate the saturation model introduced by Iancu, Itakura 
and Munier [12j | to take into account the heavy quark contribution to the proton structure function. The resulting fit 
provides a very good description of the HERA measurements of F% and the predictions for the heavy-quark structure 
functions as well as for the longitudinal structure function are in good agreement with the existing data. 
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The present work, using the dipole picture formalism, has two important features that distinguish it from previous 



studies in the literature: firstly, following the idea in 12j, we have directly use the predictions from high-energy 



QCD to parametrise the dipole-proton amplitude. The saturation properties, e.g. geometric scaling and its window 
of validity, are thus parametrised as they are predicted from perturbative QCD. Secondly, the saturation scale that 
result from our fit is not significantly reduced compared to the saturation scale obtained with light quarks only. 
This contrasts with previous studies in the literature which report a decrease of the saturation momentum. In the 
present approach saturation effects cannot be neglected. This new analysis gives an additional argument in favour of 
saturation. 

We also stress a recent analysis [26[ also including heavy-quark effects and based on the high-energy QCD properties, 
where the amplitude is parametrised in momentum space rather than in coordinate space. Compared to that study, the 
parameters for the saturation scale and geometric scaling window extension, both slightly too small in the momentum- 
space parametrisation, are more reasonable in the present approach. The analysis in momentum space was however 
carried with fixed j c and further studies are requested to check whether a higher critical slope also improves the 
model. 

A last comment concerns the relation between the parameters obtained in our fit and the predictions from the 
NLO BFKL kernel with a saturation boundary. Indeed, though a bit smaller, the exponent of the saturation scale 
remains close to the predictions from the NLO BFKL Q. In addition, we have seen that the critical slope 7 C kept 
fixed to its leading-order value in the massless case, is no longer in agreement with that value once the heavy quarks 
are taken into account. The value we obtain in that case is rather in agreement with the slope one would obtain 
from the renormalisation-group-improved BFKL kernels at NLO. This finding is another welcome outcome of our 
parametrisation. 

The last point naturally suggests to directly compare the predictions of NLO BFKL with saturation effects to 
the proton structure function. However, to rigorously address that question, one probably also needs to correctly 
introduce the running-coupling corrections to the present formalism. Though they are not expected to lead to 
significant modifications, they come with a few changes that we leave for future work. 

Finally, now that the coordinate-space dipole-proton scattering amplitude is available as directly predicted from the 
QCD saturation framework, applications to other observables are to be done. This includes the diffractive structure 
function often considered as an excellent candidate for the observation of saturation and well-described in the dipole 
picture @, [HJ [HI, H3, [28| • In addition, extending to non-zero momentum transfer (also predicted by high-energy QCD 
|29|) or parametrising the impact-parameter dependence, it can also be tested against the diffractive vector-meson 
production and DVCS cross-sections [3(| H3, H2] • Those complementary analysis are also left for future work. 
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